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The methods allowing to extract the coherent component of pion emission conditioned by the 
formation of a quasi-classical pion source in heavy ion collisions are suggested. They exploit a 
nontrivial modification of the quantum statistical and final state interaction effects on the correlation 
functions of like and unlike pions in the presence of the coherent radiation. The extraction of the 
coherent pion spectrum from tt + tv~ and tt^tt^ correlation functions and single-pion spectra is 
discussed in detail for large expanding systems produced in ultra-relativistic heavy ion collisions. 
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I. INTRODUCTION 



The hadronic observables, such as single- or multi-particle hadron spectra, play an important role in the studies 
of ultra-relativistic heavy ion collisions. However, these observables contain rather indirect information on the initial 
stage of the collision process since the particle interactions result in substantial stochastization and thermalization 
of a system during its evolution. Nevertheless, the final hadronic state can carry some residual signals of the earlier 
stages of the particle production process. A partial coherence of the produced pions is supposed to be one of the 
£C) • important examples. 

£> ' The first systematic study of coherent processes in high energy hadron-nucleus (h + A) collisions was based on 

Glauber theory [1]. In this theory, the h + A collision is considered as a process of subsequent scatterings of the 
projectile on separate nucleons of the nucleus; the projectile energies are supposed much higher than the inverse nucleus 
radius (Eh ^$> 1/R), thus allowing to consider a linear projectile trajectory inside the nucleus (eikonal approximation). 
If the scattering process occurred with almost no recoil of the nucleus nucleons, i.e. with no witnesses of the individual 
,_^ \ scatterings, then the h + A collision should be described by a coherent superposition of the elementary hadron- nucleon 
scattering amplitudes. Such a type of the collision is called coherent scattering. Since the nucleus in coherent scattering 
does not change its state, it manifests itself just as a particle with some form-factor. In the oscillator approximation, 
the nucleus form-factor can be represented by a Gaussian: exp(— q 2 i? 2 /4). The coherent processes are essential only 
for small momenta transferred from the projectile hadron to the nucleus: |q| < 1/R. Then, one can neglect the recoil 
energy and consider the nucleus as a whole during the scattering process. There is a kinematic limitation of the 
minimal longitudinal momentum transfer, |<7z| m i n ~ (M 2 — m 2 )/(2 |ph|), required to produce a particle or a group 
of particles of the invariant mass M. The vanishing of |<Zz| min with the increasing energy explains why the coherent 
processes can take place only at high enough energies. It is worth noting that the total coherent cross-section does 
not die out with the increasing energy (see., e.g., [2]). 1 

Typically, however, the transferred momenta are sufficient for substantial recoil effects and the excitation of the 
nucleus or its breakup. Then, due to a small coherence length ~ 1/ |q|, the nucleus does not participate in the 
collision as a whole and one can consider the h + A collision as an incoherent superposition of elementary hadron- 
nucleon scatterings corresponding to random phases of the amplitudes of the latter. The resulting cross-section is 
then given by the sum of the moduli squared of the amplitudes (probabilities) at each of the possible scattering points 
(unlike to coherent scattering, when the individual amplitudes are summed up first). As a result, the particles are 
produced in chaotic (incoherent) states. 

Let us come back to the production of particles (e.g., pions) in the processes of non-elastic coherent scattering at 
small transferred momenta. Since the nucleus is not excited in these processes and manifests itself as a quasi-classical 
object, one can describe particle production using the quantum field model of interaction with a classical source [3]. 
It is well known that the interaction with a classical source results in the production of bosons in coherent states [4] . 
These states minimize the uncertainty relation and, so, are the closest to classical ones. 2 This is the main physical 
link between the processes of coherent scattering and particle production in coherent states. 



1 We are grateful to V. L. Lyuboshitz for drawing our attention to this important point and for an interesting discussion. 
2 The coherent states have been introduced and studied in detail by Glauber [5]. The concept of coherent states was then 
applied to pion production in high energy processes in Refs. [6-8]. 



In heavy ion collisions at high energies, the average multiplicities are quite high, e.g., several thousands of pions can 
be produced at maximal RHIC energies. The inclusive particle spectra thus represent natural characteristics of these 
processes. A convenient way to account for the coherent properties of these processes consists in a model description of 
particle emission, rather than in detailed evaluation of the contributing amplitudes. The Gyulassy-Kauffmann- Wilson 
(GKW) model [8] is an example of such an approach. The model assumes that all pions are radiated by classical 
currents (sources) which are produced in some space-time region during the collision process. The corresponding 
density matrix is constructed by averaging over the unobservablc positions of the centers of individual sources. The 
pion spectra then effectively contain both chaotic and coherent components. In fact, the chaotic component dominates 
in case of a large emission region, while, in the opposite limit of very small space-time extent of this region, almost all 
pions are produced in the coherent state. This seems to be rather general result: if the distances between the centers 
of pion sources are smaller than the typical wave length of the quanta (the source size), the substantial overlap of the 
wave packets leads to the strong correlations (indistinguishability) between the phases in pion wave functions and, 
thus, to the coherence [9,10]. 

Recently, the coherence of multipion radiation in high energy heavy ion collisions was studied within GKW model 
in Ref. [11]. In the model, due to the longitudinal Lorentz contraction of the colliding nuclei, almost all pions 
produced with small transverse momenta p t < 1/R in central nucleus- nucleus collisions are emitted coherently, and 
their momentum spectra are determined by the system's space-time extent. Clearly, the coherence of pions can be 
destroyed by pion rescatterings. Nevertheless, the duration of hadron formation may happen to be long enough to 
allow a considerable part of the coherent pions escape from the interaction zone without rescatterings [11]. However, 
as noted in [11], one can expect a strong suppression of the GKW mechanism of coherent pion production if quark- 
gluon plasma were created: the hadronization then occurs in a thermal quark-gluon system and hadrons are produced 
in the chaotic state only. Note that clear signals of the thermalization and collective flows, observed at CERN SPS 
and RHIC energies (see, e.g., [12,13] and references therein), point to strong rescattering effects and may reflect also 
the importance of the quark-gluon degrees of freedom. 

The new physical phenomena, expected in RHIC and LHC experiments with heavy ions, are associated with the 
creation of quasi-macroscopic, very dense and hot systems. In such systems, the deconfinement phase transition 
and the restoration of the chiral symmetry arc likely to happen, possibly leading to creation of the new states of 
matter: quark-gluon plasma (QGP) and disoriented chiral condensate (DCC). In the latter case, another possibility 
for the coherent pion radiation (above the thermal background) appears. If the DCC were created at the chiral phase 

transition, a quasi-classical pion field tt c i forms the ground state of the system. The subsequent system decay is 
accompanied by a relaxation of the ground state to normal vacuum. Such a process can be described by the quantum 
field model of interaction with a classical source (see, e.g. [14]), and results in the coherent pion radiation. One of 
the general conditions of the ground state rearrangement and formation of the quasi-classical field is a large enough 
system volume [15]. Therefore, such a field could be generated in heavy ion collisions at sufficiently high energies 
provided the spontaneous chiral symmetry breaking via DCC formation takes place. The overpopulation of the (quasi) 
pion medium, making it close to the Bose-Einstein condensation point, can lead to the strengthening of the coherent 
component conditioned by the ground state (quasi-particle vacuum) decay [16]. Since the DCC appears relatively 
late (at the end of the hadronization stage), the coherent radiation could partially survive and be observed. 

The coherent emission manifests itself in a most direct way in the inclusive correlation function C(p, q) of two 
identical bosons in the region of very small |q|; p = {p\ +p^)/2, q = p\ — p2- In case of only chaotic contribution, 
the intercept of the quantum statistical (QS) Bose-Einstein part of the correlation function Cqs(p, 0) = 2 [17] while, 
in the presence of the coherent radiation, Cqs(p, 0) < 2. Generally, the coherence means strong phase correlations of 
different radiation components. In Ref. [9], a simple quantum-mechanical model of the phase-correlated one-particle 
wave packets with different radiation centers has been considered. In such a case (corresponding to indistinguishable 
correlated emitting centers), the emission amplitude A{p) averaged over the event ensemble is not equal to zero, 
(A(p)) =/= 0, and the QS correlation function intercept Cqs(p, 0) < 2. In the second quantization representation (more 
adequate for processes of multi-boson production), the analogous results take place for inclusive averages of the 
quantum field operators: (a(p)} ^ 0, Cqs(p,0) < 2, provided the radiation has a non-zero coherent state component. 
The latter represents a superposition of the states of all possible boson numbers at fixed phase relations. 

In practice, most of the correlation measurements is done with charged particles. However, charged bosons cannot 
form the usual coherent state since it obviously violates the super-selection rule. To overcome this difficulty, the 
generalized concept of charge-constrained coherent states should be used [7,8,18]. Nevertheless, the correlations of 
charged bosons are usually described with the help of ordinary (not charge-constrained) coherent states [19,20] (see, 
however, Refs. [21,22]). Our treatment of two-pion correlations takes into account the restrictions imposed by the 
super-selection rule and is based on the density matrix formalism. 

The density matrix approach gives the possibility to describe, in a natural way, the chaotic radiation (the initial 
state then corresponding to a local-equilibrium statistical operator of quasi-particle excitations) and coherent emission 



(arising due to the interaction with a classical source). This approach can easily incorporate also the squeeze-state 
component of pion radiation [23] , appearing due to the modification of the energy spectrum of quasi-pions as compared 
with that of free pions [24] . The density matrix formalism is also simply related with the Wigner function description 
of the multiparticle phase-space and its evolution governed by the relativistic transport equation [25], representing 
very useful tools with a clear classical limit. Recent development of the classical current approach to multiparticle 
production [23,19] has made it closer to the density matrix formalism; particularly, the clasical current in momentum 
space has been shown mathematically identical with the coherent-state representation of the density matrix, the latter 
called "P" or Glauber-Sudarshan representation [5], see also [26]. 

In our approach, the super-selection rule requires an averaging, in the density matrix, over all orientations of the 
quasi-classical pion source in the isospin space. As a consequence, the averaged pion field vanishes: (a(p)) — whereas, 
for identical pions, the intercept Cqs(p, 0) is still less than 2. The correlations of non-identical pions also appear to 
be sensitive to the presence of the quasi-classical source. This sensitivity arises due to properties of the generalized 
coherent states satisfying, after the averaging over all orientations of the quasi-classical source in isospin space, the 
super-selection rule for charged particles. Due to isospin symmetry of the strong-interaction Hamiltonian, there are 
unique relations for the intercepts Cq S (p, 0) of the pure QS correlation functions of two pions in various charge states 
i, j = ±, 0. For example, the coherence suppression of C determines the coherence enhancement of C^ . 

The coherence phenomena can be, however, masked by a number of effects suppressing the measured correlation 
functions. The most important among them are the decays of long-lived particles and resonances (e.g., A, K®, 77, 
rj ', ...), the single- and two-track resolution and particle contamination. In Ref. [27], the method to discriminate 
between the effects of coherent radiation and decays of long-lived resonances has been proposed. The method assumes 
the simultaneous analysis of two- and three-particle correlation functions of identical pions. The practical utilization 
of the method is however difficult due to a low statistics of near-threshold thrcc-pion combinations and the problem 
of the three-particle Coulomb interaction; also, one has to account for the super-selection rule. 3 Therefore, in the 
present work we will restrict ourselves to the consideration of two-particle correlation functions. 

In addition to QS, the correlations of particles with small relative velocities are also influenced by their final state 
interaction (FSI) . The effect of the latter on two-particle correlations is well understood and introduces no principle 
problems. It is important that the correlations in different two-pion systems are influenced by the QS, FSI and 
coherence effects in a different way. This offers a possibility to discriminate different effects suppressing the measured 
correlation functions and so to extract the coherent contribution using correlation functions of like and unlike pions 
measured at small relative momenta. 

In the paper we study the influence of the coherent pion radiation on the behavior of pion inclusive spectra and 
two-pion correlation functions and, based on it, develop the methods for the extraction of the coherent component 
above the chaotic background. Despite we associate the coherent radiation with the formation of the DCC (as the most 
probable mechanism of the coherence in ultra-relativistic A+ A collisions) , our results are rather general. Actually, 
they are based on the general properties of the coherent pion radiation: the quasi-classical nature of the coherent pion 
source and the constrains imposed by the charge super-selection rule. 

In Sec. II, we consider a general form of the density matrix of partially coherent pions, and calculate quantum 
statistical correlations of identical and nonidentical pions. 

In Sec. Ill, we set forth the density matrix formalism taking into account the decays of short-lived resonances and 
FSI of produced pions, and calculate the corresponding correlation functions. 

In Sec. IV, we discuss how to extract the coherent component of particle radiation from the two pion correlation 
functions, particularly, in the case of large expanding systems produced in ultra-relativistic A+A collisions. 

A short summary and conclusion are given in Sec. V. 

II. QUANTUM STATISTICAL CORRELATIONS OF PARTIALLY COHERENT PIONS 

It is well known that the description of the inclusive pion spectra and two-pion correlations is based on a computation 
of the following averages [8] : 



3 The latter problems are absent for neutral pions. However, sufficiently accurate measurements of neutral pion correlations 
are practically out of the present experimental possibilities. 
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where T(in;p, a) is the normalized invariant production amplitude. The summation is done over all quantum numbers 
a of other produced particles, including integration over their momenta; aj(p) and ai(jp) are respectively the creation 
and annihilation operators of asymptotically free pions i = ±, 0; the bracket (. . .) formally corresponds to the averaging 
over some density matrix |/)(/|. A special attention requires the production of particles with near-by velocities which 
can be strongly influenced particle interaction in the final state. In this Section, we concentrate mainly on quantum 
statistical correlations ignoring, for a while, the effects of resonance decays and FSI. 

Let us suppose that the density matrix p is a statistical operator describing the thermal hadronic system in a pre- 
decaying state on a hyper-surface of the thermal freeze-out 07 : t = tf (x) . After the thermal freeze-out the system 
is out of local thermal equilibrium but still can be in a pre-decaying (interacting) state. In fact, the complete decay 
(neglecting the long-time scale forces) happens at some finite asymptotic times t out < 00. Then the formal solution 
of the Heisenberg equation for the pionic annihilation (creation) operators at this post thermal freeze-out stage has 
the form 4 : 

&i,qm(P,tout) = fa,qm(P,tf) + d i (p,tf,t ut)]e~ iU ' p{tout ' tf) - (2) 

It formally corresponds to the sum of the general solution of the free (homogeneous) Heisenberg equation of motion for 
pionic field (first term), and a particular solution of the complete (inhomogeneous) Heisenberg equation with a source 
(second term). The value dj(p,t/, t out ) depends on the actual form of the source term in the Heisenberg equation. 

The decay of the system at this stage, tf < t < t ouU can be accompanied by the coherent pion radiation due to the 
modification of hadron properties in hot and dense hadronic environment or - due to some peculiarities of the phase 
transition from QGP to hadron gas, e.g., the formation of DCC. In both cases, almost non-interacting quasiparticlc 
excitations could be formed above a rearranged ground state ("condensate"). 

In the systems containing the DCC, the appearance of the quasi-classical pion field 7r c ; (corresponding to the 
density of virtual pionic excitations of the quasi-pionic vacuum) at the thermal stage is usually described in the mean 
field approximation as iri, c i(x) = ni(x) — TTi. qm (x), where the field Tri, qm (x) corresponds to the quasi-pion quantum 
excitations above the temporary vacuum background iti^ c i{x) (the order parameter). Assuming the isotopic symmetry 
of the Lagrangian like in the sigma model (see, e.g., [28]), we have Ki. c i{x) — e,%K c l{x), where e is randomly oriented 
unit vector, e 2 = 1, in the three-dimensional isospin space. Then, for each e— orientation of the quasi-pionic vacuum 
at the thermal freeze-out, the free quasi-pions 7r qm axe distributed according to the Gibbs local-equilibrium density 
matrix p e above the quasi-pionic vacuum. After the thermal freeze-out, when the decay of such a thermal system 
happens, the quasi-pion masses approach the usual free particle values and the condensate (the temporary disoriented 
vacuum) tends to relax back to the normal vacuum by emitting physical pions in coherent states - the vacuum for 
quasi-particles becomes a coherent state for free particles. The latter process is similar to particle radiation by a 
classical source. 

Then the "source" term in Eq. (2) takes on the form 

di(p,tf,tout) = di,qm(p,tf, tout) + eid co h{p,tf, tout), e =cos# , e± = — =-e ±4< ^, (3) 

where di_ qm (p,tf,t out ) and eid co h(p,tf,t ou t) are Q- an d c-value quantities respectively. While the total number of 
pions of momentum p radiated by a classical source is fixed by |d co /((p,t/,t ou t)| 7 the distribution of radiating pions 
in isospacc is determined by the orientation of the vector e; we suppose e independent of x. Wc further assume that 
the quasi-pion masses at the thermal freeze-out are near the physical mass, rrii(tf) ~ m out = m, neglecting a possible 



4 For a space-like hypersurface 07 (an example is af — i/(x) = (r 2 + xf ong ) 1 ' 2 in the Bjorken hydrodynamic model with the 
proper expansion time t), the use of the covariant Tomonaga-Schwinger formalism gives the same result with the substitution 
t -» t(x). 



mass shift which can generate squeeze-state components in particle radiation. 5 We will neglect the rescatterings at the 
post thermal freeze-out stage, i.e. put di tqm (p,tf,t out ) w 0, and approximately describe the production of coherent 
pions at this stage by the quantum field model of the interaction with a classical source [3] . Then, there is well known 
linear relationship between the annihilation (creation) operators diagonalizing the pion field Hamiltonian at the times 
tf and tout (i — ±,0): 

&i,qm(p,tout) = [&i,gm(P,tf) + eid coh {p,tf,t O ut)}e~ i0Jf,(tout ~ tf) , (4) 

where the c- value quantity d co h(p,tf,t ut) depends on a mechanism and the rate of the classical field decay. 6 

The operators a,i(p) of the asymptotic free pion field (with the origin of the time coordinate shifted to the point tf) 
are connected with the operators a, >gm (p, t) taken at the asymptotic times t ou t by the relation [30] 

diip) = Vpoe ipo{t ° ut ~ tf) 8H, qm (p, tout), Po = u> p . (5) 

Eqs. (4)and (5) allow to calculate the mean values of the asymptotic operators ctiip) and aj(p) for each e-orientation 
of the quasi-pion vacuum applying the thermal Wick theorem to the operators aj iQTO (p, tf) and a? „ m (p, tf). The 
Gaussian form of the statistical operator p e guarantees that {&i. qm (p,tf)) e = for any fixed isospin orientation e of 
the quasi-particle vacuum. Then, 



5 ij 



(6) 



(al(p2)a t (pi)) e (al(p 1 )a t (p2)) e -(al(pi)) e (al(p 2 )) e (a l (p 1 )) e (a l (p2)) 
Here 

(ot(pi)oi(p2)) e = <oJ(pi)oi(p2)) cfc + <al(Pi)) e (a 4 (P2)) e , (7) 

where the irreducible (thermal) part of the two-operator average 

(4(pi)ai(p 2 )) ch = VPwP2o(4, qm (Pl> t f)^qm(P2,tf)) e (8) 

does not depend one 7 and 

( a i(p)) e = e i d (p) = eiy/pod co h(p,tf,t out ) . (9) 

One can introduce the one-particle Wigner function [25] 

U,i{x,p) - (2^)- 3 Jd 4 q'Kl , -P)e i9 ' x (al(p + q'/2)a i (p-q'/2)) e , (10) 

satisfying the relation 

p^f e ,i(x,p)=0 (11) 

and describing the phase-space density of the non- interacting pions at t ^ t out or, in covariant formalism, at t ^ 
o~out — iowt(x); here a out is a space-time hypersurface where the interactions are "switched off" and particles can be 
considered as free. From Eq. (10), we get 



5 Squeeze-state component can arise also in a strongly inhomogeneous thermal boson system for particles with wavelengths 
larger than the system's homogeneity lengths [29]. Below we will assume the pion Compton wave-length much smaller than 
the typical system lengths of homogeneity (e.g., hydrodynamical lengths) at the thermal freeze-out hypersurface 07. 

6 It follows, from the continuity of the complete field 71"; (2:) and its derivative at t = tf that, for a fast freeze-out (t out —tf — » 0), 
the quantity d co h{p,tf,tout) is directly associated with the strength of the pion condensate. On the other hand, an adiabatically 
slow switch-off of the classical source yields d C oh(p,tf,t ou t) ~ [3]. 

Such a dependence could take place if the mass shift were non-zero and dependent on the e-orientation of the quasi-pion 
vacuum. 



(ot(pi)oi(p2)) e = / da^f e , i (x,p)e- i ^, q = pi~P2, P =(pi+p 2 )/2. (12) 



Using Eqs. (7), (9) and (12), one can split the Wigner function into the chaotic (c/i) and coherent (coh) components: 

fe,i( x 'P) = fch(x,p) + \e t \ 2 f coh (x,p). (13) 

Integrated over <J outl these components determine the operator averages (aj(pi)aj(p2)) h an d ( a l(Pi)) ( a i(P2)) 
respectively: 

(a\(pi)a t (p 2 )) ch = / dcj^e- iq - x f ch {x,p), 

(14) 
<«J(Pi)»i(P2)> e = N 2 cf (pi)d(p 2 ) = |e 4 | 2 / &r„p" e -*«-*/cofc(a;,p). 

Tout 

We suppose that the system has zero average charge and calculate the observablcs averaging over the random 
orientation of the quasi-pion vacuum in the isospin space (dfl(e) — d cos 9d(f>): 

(...) = Sp(... P ) = (47T)- 1 [dSl(e){. . .) e = (^r 1 fdn(e)Sp(. . . Pe ). (15) 



The observable pion field is related to the ensemble of events only so the corresponding complete averages of the 
asymptotically free operators vanish, for example, (a 7r +(p)) = (47r) _1 J dCl(e)(a w +(p)) = 0. The averages of these 
operators also vanish for charge-constrained coherent pion states |c), the states of a fixed electric charge and isospin - 
so called generalized coherent states [7,8,18]. This means that the density matrix p can be represented as a weighted 
sum of the projection operators \c)(c\ of these states. 

To illustrate this statement, let us consider a simple artificial case of only two sorts of oppositely charged bosons 
in one mode. Then the usual coherent states \a\), A = ±, are 

\a\) = cxp(-i \a x \ 2 ) f] f „m /2 |»A), a\\a\) = a\\a\) , 

(16) 
|n A ) = (n!)- 1 /2( a t)«|o A ) ; [a A ,4,] =8 X x>, a± = |a|e ± ^. 

These states represent superpositions of the states with different charges and so violate the super-selection rule. The 
charge-constrained coherent state |co) of charged quanta with a zero total charge may be obtained by projecting this 
state out from the charge-unconstrained two-component coherent state a+) a_) [18]: 

1 2 J °° I I 2 ™ 

|co> — / d<p\a + )\a_) = cxp(- |«| 2 ) £ L±-|„ + )|„_). (17) 

o ™ =0 

One may see that the zero charge state \cq) represents a superposition of the states with the same charges (with equal 
numbers of particles and antiparticlcs) and thus satisfies the super-selection rule. Similarly, the density matrix 

2tt 

P = 2T / d(j)\a + ) |a_)(a+| (a_| = 
o 

(18) 

l Ol |2n ^ ^ ^ ^ | a |»i+» ! i+»s+»4 I \l \/ 1/ I 

exp(-2|a| ) L L L L ( n ,l ) l/2j 2 m/2 (m m/2 (ni i ) l/2 ^i-n 2 ,n 3 -n4 \ n h + ) "2,-)(»3,+ ("4- 
ni=0n 2 =0n 3 =0n 4 =0 v ' y ' y ' y ' 

describes the mixture of the charge-constrained coherent states |c„): 

OO 

p= E i c «>< c «i' ( 19 ) 

n— — oo 

where \c n ) is the coherent state of charge "n": 



2 ~ ~ \a\ ni+n2 e 1 ^ 1 -^ 

\c n ) = exp(- \a\ ) 2^ 2^ 6 ni-n 2 , n f ,u /2r n i/2 K+)|n 2 ,-). (20) 

While, in our example, the system described by the density matrix p has not a definite charge, the average charge is 
equal to zero: 

Sp(p(a\.a+ - ala-)) = 0. (21) 

Note, that the expectation values of the annihilation operators in the corresponding coherent states are non-zero, 
(a x \a\\a\) = a\, while Sp(pa x ) = 0. 

Continuing the discussion of coherent pion production, we will assume the density matrix p e of a Gaussian-type in 
terms of the quasi-particle annihilation (creation) operators aj )gm (p,t/), related to the free particle operators according 
to Eqs. (4) and (5). Then, similar to the above example, this density matrix can be expressed through the projection 
operators on the usual charge-unconstrained coherent states of free pion field. Averaging p e over all directions of the 
isovector e according to Eq. (15), we finally get the density matrix p in a form of a weighted sum of the projection 
operators on the charge-constrained coherent states describing, in agreement with the super-selection rule, the system 
of a fixed average charge. 8 

The expressions for pion spectra in Eq. (1) thus contain the averaging over the direction of the isovector e. As a 
result, the single-pion spectra are independent of pion charges i — ±,0: 

W PW = (^)- 1 Jdn(e)Jda^f e , i (x,p) = Jda^f(x,p), 

(22) 

f(x,p) = fch{x,p) + \fcoh{x,p), 



where we have used the equality (4*7r) 1 Jdfi(e)|ej| 2 = 1/3. Note that the coherent part of the single-pion spectrum is 

d 3 N coh _ d 3 N d 3 N ch 1 /■ 1 2 

w P -j5 = ^p-^-GiP) = w p— 75 — D[p) = - / dcr^fcohKXiP) = -\d(p)\ , (23) 



where the functions G(p) and D{p) measure the coherent fraction: 

r( ) - D ^ -= d3N ^oh/d 3 p _ \$ da^fcohtxip) _ d 3 N coh /d 3 p _ § J do ^ j ' coh (x , p) 

[P) ~1 + D(p)~ d 3 N/d 3 p ~ fdWf(x,p) ' {P) ~ d 3 N ch /d 3 p ~ Jda^f ch (x,p) ■ [ ' 

The coherence influences also the quantum statistical (without FSI) correlation functions: 

C ij , } = (^y 1 Jdn(e)(4(p 1 )a}(p 2 )a t (p 1 )a 3 (p 2 )) e 

QSJK ' 1 / (4^)-i/^(e)(al(p 1 )a i ( Pl )> e ) ((4tt)-i f<ffi(e)(at(p 2 )a> 2 )> £ N ' 



Taking into account Eqs. (6)-(9), (24) and the equalities pi )2 =p± q/2, we get 

C% s (j>,q) = 1 + (9(|e 4 e,f) - 1 - 5 t] )G{ Pl )G{p 2 ) + 5 t] (cos{qx 12 ))' 

= [1 + D( Pl )] - 1 [1 + D( P2 )] _1 {l + D( Pl ) + D( P2 ) + 9(|e iej \ 2 )D{ Pl )D(p 2 ) ( 26 ) 

+8.^{cos(qx 12 ))[ h [l +V(p ll p 2 ) +V(p 2 , Pl )\Y 

where the quasi-average (cos(<7Xi 2 ))' = (cos(g(xi — x 2 )))' is defined as: 

(cosf x ))' = / rf3fJ M( a: i) rf3cr ^( a: 2)p A 'p' y /(a:i,p)/(x2,p)cos(gxi2) 
f d 3 <7^(x 1 )d 3 <r l/ (x 2 )p f fp 2 :f(x 1 ,p 1 )f(x 2 ,p 2 ) 



We do not consider here the squeeze-states of the density matrix conditioned by possible mass shift of quasi-particles. Note, 
however, that charged pions have anyway no squeeze-state components [23]. 



and similarly, with the substitution / — ► f c h, the quasi-average {cos{qx\2))' ch ; the function 

lJda^f coh (x,p)e- i<lx _ \d*{ Pl )d( P2 ) 



Note that 



P(Pi,P2)= H . ""," ^— - . ' 77,/ , D(p,p) = D(p)- (28) 

J da^f ch (x,p)e ^ x (oJ(pi)oi(p2)) cfc 



(cos( 9 xi2J) = G(pi)G(p 2 ) + h T m mm _i_ rv m cos g^iz))^ 

[l + D(pi)J[l + D(p 2 )J 

(29) 
_ 1 + D(pi,P2) + D(P2,Pi) + V( Pl ,p 2 )V(p 2 , Pl ) 



Calculating the averages 



[1 + DfcOHl + !>&*)] <«»(^)^- 

(he,| 2 ) = (47T)- 1 fdn(e)\ ei e 3 \ 2 , (30) 



<| e o| 4 ) = ^ <| e± | 4 ) = (| e+e -| 2 ) = ^, <|e e ± | 2 > = ^, (31) 

we get for the intercepts of the QS correlation functions: 

C++( P ,0) = 2 - |G 2 ( P ), C™(p,0) = 2 - iG 2 (p), 

(32) 
C&(p,0) = 1 + ±G 2 (p), G+°(p,0) = 1 - §G 2 ( P ). 

Particularly, it follows from Eqs. (32) that the decay of the quasi-pion vacuum suppresses the correlation functions 
of identical charged pions and enhances the one of non- identical charged pions, the latter effect being by a factor of 
4 smaller. For G 2 (p) = 1, the intercepts in Eqs. (32) coincide with those found in Ref. [31] in the case of a strong 
pion condensation. Our results however differ from the intercepts found in the model [21,22] of pion emission in a 
pure quantum state, - the charge-constrained coherent state. They are recovered only for large average numbers of 
coherent pions. One can then replace the canonical ensemble corresponding to the pure quantum state with a fixed 
charge, by the grand canonical one, described by the density matrix of the ensemble with a fixed average charge. For 
ultra-rclativistic A+A collisions, the inclusive description based on the grand canonical ensemble is a fairly adequate 
approach, allowing to built explicitly the density matrix for a mixture of thermal and charge-constrained coherent 
radiations and make some calculations analytically. 

One can check that the intercepts, as well as the QS correlation functions at any q, satisfy the relation [32] 

C qs + C qs = C QS + C qs- ( 33 ) 

This relation follows from the assumed isotopically unpolarized pion emission. It is valid also for the complete 
correlation functions (with FSI), except for the region of very small |q| where the correlation functions of charged 
pions are strongly affected by the isospin non-conserving Coulomb interaction. 

Note that the correlation functions, as well as their QS parts, satisfy the usual normalization condition C(p, q) — ► 1 
at large |q| provided that the coherent part of the Wigner density vanishes with the increasing \p± q/2\ faster than 
the chaotic one, i.e. G(p ± q/2) — > at large |q|. 

To get some insight in a possible behavior of the relative coherent contribution G(p), consider the situation when the 
system decays during rather short time, t out — tf — > 0, and the partial (at a fixed e) average of the pion annihilation 
operator has a simple Gaussian form: 

<a l (p)) e ^exp(-i? 2 oh p 2 ). (34) 

According to Eq. (14), the corresponding Wigner density 

fcoh(x,p) - cxp(-2i? 2 oh p 2 - x 2 /2 J R c 2 J, (35) 

so the parameter R co h determines not only the spectrum, but also the characteristic radius of the region of the 
instantaneous coherent pion emission in accordance with the minimized uncertainty relation AxAp = h/2. Let us 



assume a similar Gaussian parametrization of the chaotic component of the Wigner density in the non-relativistic 
momentum region: 

f ch (x,p) ~ cxp(-2i?2,p 2 - x 2 /2i? c \), (36) 

where Rt = (4mT) -1 / 2 measures the characteristic size of the single-pion emitter (heat de Broglic length) and 
Rch > Rt is the characteristic radius of the region of the chaotic pion emission. In the considered rare gas limit, we 
then get the correlator 

(cos(qx 12 )Y ch = cxp(-i? 2 q 2 ), (37) 

where R = (R 2 h — R 2 ^) 1 ' 2 ~ Rch represents (in the absence of the coherent contribution) the usual interferometry 
radius. The coherent fraction G(p) = D(p)/[1 + D(p)] and 

d 3 N coh /d 3 p _ U da^f coh {x,p) „ 2 21 ( no, 

D ® = ^W^F = fda^f ch ( X ,p) ~ 6XP [~ 2 ^ - R ^ P ] ' (38) 

We see that G(p) — > at large |p| on a reasonable condition R co h > Rt- 

In fact, since the quasi-classical (coherent) pion emission is conditioned by the decay of a thermal system, one may 
expect the effective radius for the coherent radiation, R co h, close to that for the thermal emission, R c h- Generally, 
in dynamical models, the effective radius varies with the momentum p and characterizes the size of the homogeneity 
region - the region of a substantial density of the pions emitted at the freeze-out time with three-momenta in the 
vicinity of p. In this case, both the coherent and chaotic radii practically coincide with the homogeneity length of the 
system. Assuming R co h ~ R c h, we have T>(pi,p 2 ) ~ V(p 7 p) = D{p) and, according to Eq. (29), 

<«»(*«*,)>' « p - D(p -1+^L- _ _ (cos( qXl2 )y ch . (39) 

One can see that (cos(gxi2))' ~ (cos(qxi2))' c h a t small |q| or, in the case of a small coherent contribution D(p) <C 1. 
Note that in the opposite case, D(p) ^> 1, a decrease of the correlation function towards unity with the increasing q 2 
is conditioned by the chaotic component (cos(qxi2))' ch starting at q 2 <~ R~ 2 ln£> 2 (0) — 4p 2 . At smaller q 2 -values, the 
behavior of the correlation function is essentially flatter due to the g-dcpcndcnce of the denominator in Eq. (39). For 
the extreme case of a pure coherent radiation, D(p) — > oo (G(p) — > 1), the function (cos(gxi2))' tends to unity at all 
q irrespective of the assumption R co h ~ Rch'- 

i , u , f d 3 <j f _ l (x 1 )d 3 <T„(x 2 )p ,1 p' y fcoh(xi,p)fcoh{x2,p)cos(qx 1 2) 

cos(gxi 2 )) -► jr-^ — -——= — ji — — — — r = 1. (40) 

J d- i cr^(xi)d i a u (x2)PiP2Jcoh{Xi,Pi)J co h(X2,P2) 

The last equality in Eq. (40) follows from the definition (14) of the coherent Wigner function, both the nominator 
and denominator in Eq. (40) being equal to |<i(pi)<i(p2)| 2 - Experimentally, the approach to such an extreme regime 
can display itself as a tendency of the intercepts of the QS correlation functions to the values defined by Eqs. (32) 
at G{p) — > 1, and - as a flattering of the QS correlation functions within a growing q-interval. The latter mimics a 
decrease of interferometry radii; of course, it does not mean that the real size of the system tends to zero. 

The effect of coherent radiation on pion spectra and tt + tt + and 7r + 7r~ correlation functions is demonstrated in 
Figs. 1-3 for different ratios D to t = D(0)(Rt / R co h) 3 of the total numbers of coherent and chaotic pions. The plots 
correspond to simple Gaussian Wigner functions (35), (36) with Rt = (4mT) _1 / 2 s=a 0.72 fm (T = 0.135 GcV) and 
Rcoh = Rch = 5 fm. Under the assumption of a common source of coherent and chaotic pions in ultra-relativistic 
heavy ion collisions, characterized by a typical radius R ~ 5—10 fm, the coherent component in the spectra is 
concentrated in rather small momentum region of a characteristic width (2i?) _1 ^20 — 10 MeV/c (see Fig. 1). 
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FIG. 1. The single-pion momentum spectra d 3 N/d 3 p calculated for different ratios Dtot of the total numbers of coherent and 
chaotic pions, assuming the Gaussian parametrization of the Wigner densities in Eqs. (35), (36) with Rt = (4mT) -1 ' 2 fa 0.72 
fm (T = 0.135 GeV) and R co h = R c h = 5 fm. The solid, dotted, dash-dotted and dashed curves correspond to D to t — 0, 0.01, 
0.1, and 1 respectively. The overall normalization is arbitrary. 
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FIG. 2. The pure QS correlation functions Cqs(p, q) calculated for tv + tt + pairs at p = GeV/c on the same conditions as 
in Fig. 1. 
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FIG. 3. The pure QS correlation functions Cqs(p, q) calculated for tt + -k pairs at p = GeV/c on the same conditions as 



in Fig. 1. 



III. CORRELATION FUNCTIONS AFFECTED BY FINAL STATE INTERACTION AND COHERENCE 



In ultrarelativistic A+A collisions, free hadrons appear mainly at the late stage of the evolution after the system 
expands and reaches the thermal freeze-out. After the hydrodynamic tube decays and produces final particles and 
resonances, particles still appear from resonance decays. Thus, more than half of pions produced in high energy 
heavy ion collisions is of the resonance origin. As a consequence, the pion spectra and correlations are influenced by 
resonance production and decay spectra, as well as - by resonance lifetimes. Particularly, the pions from the decays 
of long-lived resonances do not contribute to QS and FSI correlations and thus suppress the correlation function 
C tJ (p, q) ; we will consider this suppression in next Section. 

However, even after the thermal (hydrodynamic) system and short-lived resonances decay, the particles in near-by 
phase space points continue to interact. Due to a large effective emission volume in heavy ion collisions, the particle 
interaction in the final state is usually dominated by the long-range Coulomb forces. To calculate the FSI effect on 
two-particle spectra, we will assume sufficiently small phase space density of the produced particles and use the FSI 
theory in the two-body approximation [8,33,34] for pions, neglecting the FSI of resonances. 

The single-pion spectrum in Eq. (1) then remains unchanged while the two-pion one (for pairs containing no pions 
from long-lived sources) takes the form 



^Pl^P2 



d 3 pid 3 p 2 



= /d 4 A :i d 4 fc 2 d 4 fcW 4 fc^ a l(fc 1 ) a ](fc 2 ) aj (fci)a J (^))^7^'(fc 1 ,fe)$^'*(fc / 1 ,^ 



(41) 



12 



where the non-symmetrized Bethe-Salpeter amplitude $p7p2 J (^l; k 2 ) = ^pipl''*(ki, k 2 ) in four-momentum representa- 
tion is expressed through the propagators of particles i and j and their scattering amplitude T%j analytically continued 
to the unphysical region [33,34] : 9 



K 1 ) P2 3 (h,k 2 ) = S(k 1 - Pl )S (k 2 - P2 ) + 5 (fci+fc2-Pi-P2)-^- (fc? _ roa 3 _. 0)( ^_ ro2 _. 0) . (42) 

The averaging in Eq. (41) is performed with the help of the statistical operator p without FSI: (...) = Sp(. . .p). 
Introducing the Bethe-Salpeter amplitudes ^ Pl p 2 J (x\, x 2 ) in space-time representation: 

*fc£ J '(fci,fe) - (2tt)- 8 J d 4 x 1 d 4 x 2 e ik ^ +ik2X2 ¥-^(x 1 ,X2), (43) 

one can rewrite Eq. (41) as 

d®N- ■ f 

^^' cPpKPpa " 7 ^i^^^i^^^i.^^i'^)*^'^!,^)*^'* K,4), (44) 

where the space-time density matrix p lJ is just the Fourier transform of the four-operator average in Eq. (41) 



10 



p ij (x 1 ,x 2 ;x' 1 ,x' 2 ) = (2tt)- 16 / d 4 fcid 4 fc 2 d 4 fc;d 4 fc^e 4fclKl+4fc2X2 e- tfc i<-^4( a t( fcl ) a t( fc2 ) ai ( fc ^ aj ( fc ^^ (45) 



Separating the phase factor due to free motion of the two-particle c.m.s.: 

(46) 
Xi 2 = \{xi +x 2 ), xi 2 =xi-x 2 , P = 2p = pi+p 2 

and integrating over the pair c.m.s. four-coordinates X\ 2 and X' 12 m Eq. (44), one can express the two-particle 
spectrum through the reduced space-time density matrix pp(xi 2 ;x' 12 ), the latter depending on the pair total four- 
momentum P and the relative four-coordinates of the emission points only: 



v Pl^P2 



' f ' A;; = I d 4 X 12 d 4 x[ 2 4(x 12] x[ 2 )^ q -^(x 12 )^-^%x' 12 ), (47) 



d 3 pid 3 p 2 



4(x 12 ;x' 12 ) = (2tt)- 8 / d 4 fcid 4 fcie < ( fcl - , ') x "e- < ( fc i- p > a! '«(ot(fci)ot(P-fci)oi(fci)o i (P-fci)). (48) 



Note that in the two-particle c.m.s., where P = {mi 2 , 0, 0, 0}, q — {0, 2k*}, x\ 2 = {£*, r*}, the reduced Bethe-Salpeter 

amplitude ipq (^12) = V4 (^12) at t* =t\ — t\ = coincides with a stationary solution V>-k* (r*) of the scattering 
problem having at large distances r* the asymptotic form of a superposition of of the plane and outgoing spherical 
waves (the minus sign of the vector k* corresponds to the reverse in time direction of the emission process). This 
amplitude can be substituted by this solution (equal time approximation) on condition [34] \t*\ <C mr* 2 which is 
usually satisfied for particle production in heavy ion collisions. 

Since the resonances have finite lifetimes, their decay products arc created in an essentially four-dimensional space- 
time region. At the post thermal freeze-out stage, the resonances are usually described by semiclassical techniques; 
they are considered as unstable particles moving along classical trajectories and decaying according to the exponential 
law [35] (sec, however, [33,36,37]). This approximation neglects a small correlation effect in pairs of unlike pions 
appearing due to QS correlations of identical resonances. The resonances are supposed to be described according 



9 It is important that the relation between the production amplitude and the operator product average, as given in Eq. (1), 
is valid also off mass shell. 

10 For identical particles, it differs from the space-time density matrix of ref. [33], where the effect of QS enters through the 
symmetrization of the Bethe-Salpeter amplitudes while, here - through the Wigner decomposition of the four-operator average 
in Eq. (52) below. 
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to the Gibbs density matrix prior to the thermal freeze-out; this guarantees the chaoticity of the decay pions. 11 
Therefore, the pions from resonance decays do not destroy the structure of the decomposition of the operator averages 
in Eqs. (6) and (7) into irreducible parts based on the thermal Wick theorem. 

After the production, the pions in near-by phase space points, chaotic as well as coherent ones, undergo a long-time 
scale interaction in the final state. According to Eqs. (44) or (47), the intensity of FSI interaction is conditioned by 
the two-particle Bcthc-Salpeter amplitudes ^ PlP2 (xi,X2) or tp q {x\_2) and the corresponding two-particle space-time 
density matrices p{x\, X2',x' 1 ,x 2 ) or pp(xi2',x' 12 ). Clearly, in the case of absent FSI, the two-pion spectrum merely 
reduces to the Fourier transform of the space-time density matrix. It can be represented as an integral over the mean 
four-coordinates x — (x + x')/2 of a combination of bilinear products of single-particle chaotic and coherent emission 
functions g c h(x,p) and g C oh{x,p), respectively defined in Eqs. (50) and (51) below. 

The emission function g(x,p) is closely related with the Wigner phase space density f(x,p) at asymptotic times 
t Js tout- Let us denote by x = {£, x — (p/po)(t — i)} the space-time point, starting from which a free particle moving 
with velocity p/po reaches a point x; the portion of such particles is g(x,p). Collecting all the contributions (starting 
in our case from the thermal freeze-out time tf), we have 

Pof(x,p) = j d 4 x8 3 (x- x+ (p/p Q )(t - t))g(x,p), (49) 

where g(x,p) = poS(i — tf)f(x,p) + s(x,p) and s(x,p) ~ 9(t — t)6(i — tf) is the density of pion emission at the post- 
thermal stage, t > tf. Therefore we can rewrite the irreducible (thermal) part of the two-operator average through 
the chaotic emission function as: 

(al(pi)ai(p2)) ch = J da t ,p^e- iqx f ch (x,p)^Jd 4 xe~ l i s g ch (x,p), 

(50) 
p = P/2 = (p!+p 2 )/2, q=pi-p 2 , 

where we have used the equality qx — qx following from the relation qp = qopo — qp = 0. Similarly, for the coherent 
component of the two-operator average at fixed e, we get 

<«l(Pi)) e (« 4 (P2)) e = |e 4 | 2 ci*(Pi)MP2) = 

(51) 
N 2 / d<7^e-^f coh {x,p) = \ ei \ 2 fd 4 xe-^g coh (x,p). 

"out 

The results of Section II can thus be rewritten in terms of the emission functions in accordance with a formal 
substitution J da^p^ 1 f(x,p) — > / d 4 xg(x,p). 

lout 

To express the four-operator average in Eq. (48) through the emission functions, we can exploit the decomposition 
similar to that in Eq. (6): 



<% 



<ot(fc!)ot(P - fci)oi(fci)«^ - fci)) e = <ot(fci)oi(fci)).(o5( P - fc i) fl ^ P - fc i))e+ 
"(aJ(fci) ai (P - fci)) e (oJ(P - fci)0i(fci)) e - (oJ(fc 1 )) e (ot(P - fc 1 )) e <a 4 (fci)) e («,(P - K)) t 
Using Eqs. (50) and (51) for the two-operator averages in Eq. (52), we get: 
(oJ(fci)ot(P - k^ik'Ja^P - k[)) e = fd i x 1 d 4 x 2 x 

| e _ i(fel _ fei )-X 12ffe; . ^ 1 (^ + ^ g e . ^ p _ 1 (ki + J^ + 

5 . je -i(k 1+ k[-p)-x 12 [ 5ej . ( Sl>p+ i( &1 _ fc j)) 9ei {x 2 ,p-\{kx - K)) 

-\ei\ 4 g C oh (xi,p+ \{ki - k[)) g coh (x 2 ,p~ \{k x - A;i))]}> 



(52) 



(53) 



11 Note that the chaotisation of decay pions partially happens irrespective of the form of the density matrix if pions were 
emitted by a large number of many different sorts of resonances. 
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where x± 2 — x\ — x 2 and 

g e>i (x,k) = g ch (x,k) + \e i \ 2 g coh (x 1 k). (54) 

After the averaging over the orientation of the isospin vector e, we get 

(aJ(fci)at(P-fc 1 )a i (fci)%(P-fci)> = /rf 4 Sid 4 S 2 - 

jg-i^-ki)-*" [ ff ( Slj i (fcl + fc / } ) 5 ( i2) P _ i (fcl + fc / } ) 

+ «l^ e j| 2 ) - \)9coh (xi, i(fci + k[)) g coh (x 2 ,P - |(fci + *£))]+ (55) 

^. e - i(fel + fe ;-p).x 12 [ 5 ^ UP + |(fci - fci)) 5 (x 2 ,p - §(fci - fc i)) 

-\dcoh (xi,p+ |(fci - fci))_g co/l (x 2 ,p- |(fci - fci))] |, 

where 

g(x,k)=g ch (x,k) + -g coh (x,k). (56) 

Inserting expression (55) for the four-operator average into Eq. (48) and, integrating in the first and second term over 
(fci — k[) and (fci + k[ — P) respectively, one can rewrite the reduced space-time density matrix as: 

Pp (xi2', x' 12 ) = (2tt) -4 / d 4 xid 4 x 2 d 4 K- 

|gJK-(x 12 -x' 12 )^4^1 ( Xl2 + ^ _ Si2 ) ^ (x 1)P + «) 5 (x 2 ,P - ft) 

+ {(\ e i e 3?) - \)9coh {xi,p + n)g coh (x 2 ,p- «)] + (57) 

V i "- (a,12+x " ) $ 4 (!(a;i2 - *i 2 ) - S12) [ff(*i,P + «)ff(S2,p- «) 
-^gcoh(xi,p + n)g coh (x 2 ,p- «)] J. 
According to Eq. (47) and using the equality ip q {—Xn) = ip-q{xi2), the two-pion spectrum then becomes: 

^Pi^P2 rf W^ 2 = (2Ti)- i j-d i x l d i x 2 d i nd 4 ee^- 
{[g (xi,p + k) g (x 2 ,p - k) + ((|eie,-| 2 ) - \)g C oh {xi,P + n) g coh (x 2l p - n)]ip q ~ )l3 {x 12 + \e)ij) q ~ )l3 * {x 12 - |e) 

+Sij[g (xi,p + k) g (X2,p - k) - \g co h {x\,P + k) g coh (x 2 ,p- k)]w (£12 + ^)wS q (x 12 - \e)\ 

= (2tt)- 4 / d i x 1 d A x 2 d A Kd A e e 1 ^- (58) 

\[g c h(xi,p + K)g c h(x 2 ,p- k) + l(g c h{xi,p + K)g co h{x 2 ,p- n) + g co h {x~i,P + n) g ch (x 2 ,p- «))]• 

[4~ )ij (xi2 + h)4~ ]ii *{^2 h) + s^4~ )l1 (X12 + h)4-- q )ij *(x 12 |C)] 

+ (\eiej\ 2 )g coh (xi , p + k) g C oh (x 2 ,p- «) ipq ' v (xi 2 + §e)v4 *(#12 - 5 e )j- 
If the FSI were absent, i.e. ij)q {x\ 2 ) = exp(— ig • S 12 /2), one would get 
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(59) 



"pi^ /pA = / d4 xid 4 x 2 J 5 (xi,pi) g (x 2 ,_p 2 ) + ((leiejl 2 ) - ±).g co h (xi,pi) g co h (x 2 ,p 2 ) 

+% [5 (»i . P) 5 (»2 , P) - gScofc (Si , p) g co h {X2 , p)] cos(qxi 2 ) J 

= J d 4 Xid 4 X2 |g(5i,pi)g(S 2 ,p 2 ) + ((le^ejl 2 ) - |(1 + Sij))g co h(xi,Pi) g C oh(x2,P2) 

+5 i3 g (xi,p) g (X2,p) cos(qx 12 )\ 
= J d 4 xxd 4 x 2 |g c / l (5i,pi)g ch (S 2 ,p 2 ) + i(g ch (#i,Pi) # C o/i (x 2 ,p 2 ) + g coh (xi,pi) g c/l (x 2 ,p 2 )) 

+(|e i e : ,-| 2 )g co / l (xi,pi)g co / l (x 2 ,p 2 ) + 5»j [fif c /j (xi , p) g c h {x2,p) + \gch {x\,p)g co h (x~2,p)] cos(qx 12 )\ 

and recover Eqs. (26) for the pure QS correlation functions. 

In the case of absent coherent emission, i.e. d = g co h = 0, and on the usual assumption {R^ <C Rch) °f sufficiently 
smooth four-momentum dependence of the chaotic emission function g c h(x,p) as compared with a sharp g-dependence 
of the QS and FSI correlations (determined by the inverse characteristic distance between the emission points) , the 
chaotic emission functions in Eq. (58) can be taken out of the integral over k at small values of k, this integral thus 
being close to 5 4 {e). Choosing the momentum arguments in (^-functions in accordance with Eq. (59) for the case of 
absent FSI, we get for the two-pion spectrum and the correlation function: 



« "a ~ r 74- j4- 
^pi^ Ww/p, ~Jd xid x 2 - 

\9ch{xi,Pi) g c h(x2,P2) \ipq ■ 7 (Si 2 )| 2 + Sijg ch (xi, p)g ch (x 2l p) i>q (x 12) ip Lq V * (x 12) j , 



(60) 



^«<l^ )4J (Si 2 )| 2 ) c , + %(^^"(x 12 )^*(x 12 ))^, (61) 

where the average (A) c h and quasi-average (A)' ch arc defined as: 

/ .v _ j d 4 x 1 d 4 x 2 Ag ch (x 1 ,p 1 )g ch (x 2 ,P2) , ■, 

fch j d 4 x 1 g ch (xi,p 1 ) j d 4 X2g c h{x2-,P2Y 

,j\i _ j d 4 xid 4 x 2 Agch(xi,p)g c h(x2,p) . . 

' c/l J d 4 xig ch (xi,pi) J d 4 x 2 gch(x2,P2)' 

In the case of a nonzero coherent contribution, the e/2- and Xi 2 -dispersions in the pure coherent term in Eq. (58) 
are the same (2i? 2 o/l ), contrary to usually negligible e/2-dispcrsion in the pure chaotic term: 2R\ <C 2_R 2 h . As for 
the mixed term, the e/2-dispcrsion would be negligible if only the characteristic size R co h of the coherent source were 
sufficiently small; with the increasing R co h, this dispersion may become important - for R co h ~ R c h it amounts to 
about half of the x\ 2 dispersion. Therefore, the e dependence of the Bcthc Salpctcr amplitudes should be generally 
retained in these terms. The important exception is the case of practical interest in heavy ion collisions, when the two 
charged pions are created in their c.m.s. at a distance much larger than the corresponding s-wave scattering length 
(of a fraction of fm) and much smaller than their Bohr radius (of 387.5 fm). The two-pion FSI interaction at small q 
is then dominated by the Coulomb FSI and depends only weakly on the space time separation of the emission points. 
In this case, 

Cl « (K-^12)! 2 ) + %<^- )y '(S 12 )^ )y '*(x 12 ))' + (9(| e , ej | 2 ) - 1 - 8 t3 )G( Pl )G(p2){\4-^(x 1 2)\ 2 ) coh , (64) 

where the averages are defined as in Eqs. (62) and (63) with the substitutions g c h —* g or g c h — > g C oh and, the relative 
coherent contribution G(p) - in Eq. (24) with a formal substitution / da M p M / (x , p) — > J d 4 xg (x,p). 
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IV. EXTRACTING COHERENT COMPONENT OF PARTICLE RADIATION 

Up to now, we have ignored the contributions d 3 N^ /d 3 p arising in the pion spectra from the decays of long-lived 
(I) sources such as 77-, ry'-mesons, and also the unregistered kaons and hyperons. The pions from these sources possess 
no observable FSI (due to very large relative distance of the emission points) as well as no noticeable interference effect, 
because the corresponding correlation width is much smaller than the relative momentum resolution q m { n of a detector. 
Therefore the measured correlation functions, defined in Eq. (1), can be expressed through the correlation functions 
C y (p, q) (discussed in previous Section) of all pion pairs Tr 1 !: except for those containing pions from long-lived sources 
as follows [38]: 



C^{p,q) = n ij (puP2)/n i (p 1 )n j (p 2 ) = ^{p)C^{p,q) + 1 - A«(p), 



(65) 



where the suppression parameter A lJ (j)) measures the fraction of pion pairs containing no pions from long-lived 



sources 



,.12 



A«(p) 



d 3 N % /d 3 p 



d 3 NPfd 3 p\ 
d 3 Nj/d 3 p J K L 



(66) 



In the (artificial) case of absent FSI effect, the correlation function C^(p,q) = C l A s {p,q), and the averaging in 
(cos(gxi2))' in the QS correlation functions in Eqs. (26) should be applied only to the pion pairs containing no pions 
from long-lived sources. Then, assuming sufficiently good detector resolution, q min <C R~ , we can determine the 
intercepts C lJ (p, 0) calculating the correlation functions at \q\ ~ q m in'- 



C ll { Pl ftnin) = 1 + A ij (p) Sij + (9<|e 4 e,f) - 1 - 5 tl )G 2 (p) 



(67) 



The intercepts are lower than 2 for any system of identical pions and they are higher (lower) than 1 for ir + Tr~ (it tt ) 
systems. 

Since the suppression parameters A(p) are generally different for different pion pairs, e.g., due to different contri- 
butions of hyperon decays, it is impossible, using only apparent intercepts in Eq. (67), to separate the contributions 
of the coherent and long-lived sources, unless there is known a ratio of the suppression parameters A(p) for identical 
and non-identical pions: A" (p) / A 11 (p) . Then, for example, from the intercepts of the tt + tt + and 7r + 7r _ correlation 
functions, one obtains the coherent fraction squared: 



,2 f ^ A++( P ) 



G\p) 



A+-(p) 



4 A++(p) l C++(p,g min )-l 
5A+-(p) 5C+-(p,g min )-l 



(68) 



In fact, the knowledge of the ratio A" (p) j A 1 - 7 (p) is not of principle importance for the extraction of the coherent 
fraction G(p). Besides the intercepts, one can exploit also the q dependence of Cqs(p, q) in sufficiently wide interval to 
follow Eqs. (26), and perform simultaneous or separate fits of the correlation functions C*- 7 , suitably parameterizing 
the correlator (cos(gxi2)) and the function G(p ± q/2). For example, one can use the usual Gaussian correlator 
parameterization 



(cos(qx 12 )Y ch ~ exp(-qlR 2 x - q^ - q 2 z R 2 z ) 



(69) 



in the longitudinally comoving system (LCMS) in which the pion pair is emitted transverse to the collision axis 
(Pl = 0). The components of the vector q are chosen parallel to the collision axis (z=Longitudinal), parallel to the 
vector p t (x=Outward) and perpendicular to the production plane (z,x) of the pair (y=Sideward). Assuming the 
same radii also for the coherent emission region, and a transverse thermal law exp(—m t /T) for the chaotic radiation 
with the temperature T (m t is the pion transverse mass), we can parameterize the coherent fraction G(p) similar to 
Eq. (38) for the non-relativistic case with [16] 



D{p) ~ D(0) exp 



-2(plRl 



Pytty 



v\R\) 



mi 
~T 



(70) 



12 One can include in N\ and the corresponding suppression parameters A lJ the contribution of misidentified particles which 
also introduce practically no correlation. 
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and use Eq. (39) to calculate (cos(gxi 2 ))'. 

The presence of the FSI effect introduces the additional g-dependence of the correlation functions and thus improves, 
in principle, the accuracy of the extraction of the coherent contribution G(p). Consider, for example, only effect of the 
Coulomb FSI and assume that the emission functions, g c h and g C oh, are localized in the regions of characteristic sizes 
much smaller than the two-pion Bohr radius \a\ = 387.5 fm so that the modulus of the non-symmetrized Coulomb 
wave function can be substituted by its value at zero separation. As a result the Coulomb effect factorizes in a form 
of so called Gamow or Coulomb factor A c (ak*) = \ip'?' ul {Q)\ (see, e.g., [8]): 

C(p, q) = A c (ak*)C QS (p, q), A c (x) = (2n/x)/[exp(2ir/x) - 1], (71) 

where k* = |q*|/2 is momentum of one of the two pions in their c.m.s. For the correlation functions of like (a = \a\) 
and unlike (a = — \a\) charged pions, we get 



C ±± {p,q)=A ±± {p)A c {\a\k*) 
C + -(p,q)=A + -(p)A c (-\aW 



1 + {cos(qx 12 )Y - \g{ P + q/2)G(p - q/2) 
5 



l -G(p + q/2)G(p-q/2) 



A+"(p)]. (72) 



Similar to the case of absent FSI, we can again use the parameterizations (69), (70) and the relation (39), and 
fit, simultaneously or separately, the correlation functions of like and unlike charged pions according to Eqs. (72). 
Moreover, the known g-dependence of the Gamow factors allows to separate the coherent fraction G(p) from the 
suppression parameter A(p) in a model independent way, without exploiting the q-dependence of (cos(qxi2)) c h and 
G(p± q/2). Indeed, one can perform the fits according to Eqs. (72) in an interval of q m [ n < \q\ <C R' 1 guaranteeing 
(cos(qx)}' i=a 1 and G{p\.2) ~ G(p). The q-dependence of the correlation functions is then uniquely determined by the 
known functions A c (|o|fc*) and yl c ( — |a|/c*), and the three fitted parameters: G(p), A ±=t (p) and A" 1 (p). Of course, 
such an analysis requires very good detector resolution and its good understanding. 

Note that Eqs. (72) are not applicable for very small (~ 1 fm) as well as for large sources. In the former case one 
has to account for the strong FSI, in the latter - for the finite-size Coulomb effects. For ultra-relativistic heavy ion 
collisions, the strong FSI effect on two-pion correlation functions is negligible for like charge pions and small (a few 
percent) for unlike pions. The Coulomb finite-size effects can be approximately taken into account, substituting the 
Gamow factor A c (ak*) in Eqs. (17) by the finite size Coulomb factor A c (ak*, (r*)/a) [39]. The latter represents a 
simple function of the arguments ak* and (r*)/a, where (r*) is the mean distance of the pion emission points in the 
pair c.m.s., corresponding to a given momentum p. Particularly, A c = j4 c (afc*)[l + 2(r*)/a] at k* <^ l/(r*). 

The dependence of the Coulomb factor on the unknown parameter (r*) somewhat complicates the model- 
independent method for the extraction of coherent component G(p) exploiting only the correlation functions in the 
region of very small relative momenta. Now, the simultaneous analysis of the correlation functions of like and unlike 
charged pions is required because their separate analysis yields the coherent contribution G{p) up to a correction 
(r*)/a only. As for the method based on a fit in a wide |q|-interval, the quantity (r*) being a unique function of 
the parameters characterizing the emission density, actually represents no new free parameter. Particularly, for a 
universal anisotropic Gaussian r*-distribution of the chaotic and coherent emission functions, the quantity (r*) can 
be expressed analytically through the Gaussian interferometry radii R y , R z and R* = ^jR x (M and M t are the 
two-pion effective and transverse masses respectively) in the case of practical interest, when R%> R y ~ R z [39]. 

In practice, however, the Gaussian parametrization of the relative distances between the emission points may happen 
to be insufficient. Particularly, it can lead to apparent inconsistencies in the treatment of QS and FSI effects because 
the latter is more sensitive to the tail of the distribution of the relative distances. If, for example, the r*-distribution 
were represented by a sum of two Gaussians with essentially different mean squared radii, the r*-"tail", determined by 
the larger Gaussian radius, would influence the observed correlation functions in different ways. For identical pions, 
the "tail" results in an additional rather narrow peak in the QS correlation function; however, this "tail" would show 
up only as a suppression of the correlation function if the peak were concentrated at q < q-min or if one measured 
a given projection of the correlation function (e.g., in q S id e direction) fixing others {qi ong and q ou t) in the interval 
exceeding the width of the narrow peak. At the same time, the r*-"tail" would influence Coulomb correlations at 
small q > q m [ n since the long-distance nature of Coulomb forces leads to the observable effect conditioned by the 
"tail" up to r* ~ \a\. In such a situation, one can no more rely on the equality between (r*}Qs, determined by 
the interferometry radii, and the characteristic size (r*)c determining the Coulomb FSI effect. Generally, one has 
to introduce also different suppression parameters Aqs < Ac corresponding to (t*)qs < {r*)c- Eqs. (72) for the 
correlation functions of like and unlike charged pions, with the substitution of the Gamow factor A c (ak*) by the 
finite-size Coulomb factor A c (ak* , (r*)/a) [39], are then modified to the form: 
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C ±± (p,q)=A±±(p)A c (\a\k*,(r%±/\a 



QSUV- fl cV|«'|'* , v ,q S 

±±/„u r\„\i„* /„*\±±/i..u in i 



(cos(gx))'-^G(p + g/2)G(p-g/2) 
5 



A^CpJ^dolfc*, (r*)^/|a|) + [1 - A^{p)}, 

(73) 

C+-(p, g) = A+ s (p)A c (-\a\k*,-(r% s /\a\)lG(p + q/2)G(p - g/2) + 



•5 



A+-(p)A c (-|a|fe* , -<r*)+-/la|) + [1 - A+"(p)]. 

To simplify the analysis, one can neglect a small difference between the suppression parameters Aqs and Ac due to 
the tail of the redistribution and also neglect a presumably small difference between (r*) ±± and (r*)" 1 . 

Note, that at SPS and RH1C energies the effect of strong FSI on 7t + tt~ correlations is still quite noticeable and, 
when neglected, it can lead to a suppression of a fitted (r*)+~ by ~ 50%. Also, due to a substantial inaccuracy of 
the Coulomb factor A c (ak* , (r*)/a) near the tailing point k* ~ l/(r*), the parameters (r*) ++ and (r*) H can be 
respectively overestimated and underestimated if the fitted region were not sufficiently wide. Further, in the case 
of different chaotic and coherent emission volumes, one has to use finite size Coulomb factors with different (r*) in 
the chaotic, coherent and mixed terms. All these problems can be overcome exploiting the exact formulae for the 
two-pion wave functions (in the equal time approximation) and calculating the correlation functions according to the 
approximate Eq. (64). To control the systematic errors due to the smoothness assumption in Eq. (64), one can give 
up this assumption (at least in the pure coherent term) and check the results using instead the general expression for 
the two-pion spectrum in Eq. (58). 

After the extraction of the fractions G(p) and A ++ (p) or A (p), one can obtain the coherent part of the measured 
single-pion spectra ui p d 3 N±/d 3 p. Using Eq. (66), and substituting d 3 N/d 3 p — * (d 3 N±/d 3 p— d 3 N± /d 3 p) in Eq. (24), 
one gets: 



u;- 



"If s | |d(p)|2 = -p^G^TA^). (74) 



p 



The coherent part of the observed spectra is thus directly connected with the intensity \d(j>)\ of the quasi-classical 
source of coherent pions. 

V. CONCLUSIONS 

Using the density matrix formalism, satisfying the requirements of the isospin symmetry and the super-selection 
rule for generalized coherent states, and accounting for the final state interaction in the two-body approximation, we 
have developed methods allowing one to study the coherent component of pion radiation which, in heavy ion collisions, 
is likely conditioned by formation of a quasi-classical pion source. 

These methods are based on a nontrivial modification of the effects of quantum statistics and final state interaction 
on two-pion correlation functions (including those of non-identical pions) in the presence of a coherent pion radiation 
generated by the decay of the quasipionic ground state ("condensate"), ft has been shown that the combined analysis 
of the correlation functions of like and unlike pions gives the possibility to discriminate between the suppression of the 
like-pion correlation functions conditioned by the coherent pion component and that due to the decays of long lived 
sources. 

The methods allowing to extract the coherent pion component from 7r + 7r~ and it^tt correlation functions and 
single-pion spectra have been discussed in detail for large expanding systems produced in ultra-relativistic heavy 
ion collisions. For such systems, the two-pion final state interaction is dominated by the Coulomb one and plays 
an important role in this analysis, allowing one to determine the coherent fraction using a suitable model for the 
coherent and chaotic emission functions and fitting the corresponding correlation functions. For rough estimations 
the procedure can be substantially simplified accounting for the finite-size Coulomb effects in an approximate analytic 
form [39]. 

Finally, the coherent fractions extracted from the correlation analysis, combined with the single-pion spectra, 
give us the possibility to determine the spectrum of the coherent pion radiation above the thermal background and, 
therefore, to estimate the quasipionic condensate at the pre-decaying stage of the matter evolution and discriminate 
between possible mechanisms of coherent production in ultra-relativistic A+A collisions. 
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